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ABSTRACT 

Aims. We have carried out an accurate investigation of the Ca isotopic composition and stratification in the atmospheres of 23 
magnetic chemically peculiar (Ap) stars of different temperature and magnetic field strength. 

Methods. With the UVES spectrograph at the 8 m ESO VUT, we have obtained high-resolution spectra of Ap stars in the wavelength 
range 3000-10000 A. Using a detailed spectrum synthesis calculations, we have reproduced a variety of Ca lines in the optical and 
ultraviolet spectral regions, inferring the overall vertical distribution of Ca abundance, then we have deduced the relative isotopic 
composition and its dependence on height using the profile of the the IR-triplet Ca ii line at A 8498 A. 

Results. In 22 out of 23 studied stars, we found that Ca is strongly stratified, being usually overabundant by 1.0-1.5 dex below 
logrjooo ~ -1, and strongly depleted above logrjooo = -1.5. The IR-triplet Can line at A 8498 A reveals a significant contribution 
of the heavy isotopes "^^Ca and ''**Ca, which represent less than 1 % of the terrestrial Ca isotopic mixture. We confirm our previous 
finding that the presence of heavy Ca isotopes is generally anticorrelated with the magnetic field strength. Moreover, we discover that 
in Ap stars with relatively small surface magnetic fields (< 4—5 kG), the light isotope ''"Ca is concentrated close to the photosphere, 
while the heavy isotopes are dominant in the outer atmospheric layers. This vertical isotopic separation, observed for the first time for 
any metal in a stellar atmosphere, disappears in stars with magnetic field strength above 6-7 kG. 

Conclusions. We suggest that the overall Ca stratification and depth-dependent isotopic anomaly observed in Ap stars may be 
attributed to a combined action of the radiatively-driven diffusion and light-induced drift. 
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1. Introduction 

■ After the pioneering work by Michaud (119701) . atomic diffusion 
' in stellar envelopes and atmospheres has been recognized as the 

main process responsible for the atmospheric abundance anoma- 
' lies in the peculiar stars of the upper main sequence. Due to 
, their unique characteristics, such as an extremely slow rotation, 

■ strong, global magnetic fields, and the absence of convective 
mixing, magnetic chemically peculiar (Ap and Bp) stars exhibit 
the most clear manifestation of diffusion effects and thus repre- 
sent privileged laboratories for investigation of chemical trans- 
port processes and magnetohydrodynamics. 

Detailed diffusion calculations performed for a set of chem- 
ical elements in the atmospheres of magnetic peculiar stars pre- 
dict separation of chemical elements over the stellar surface 
and with height in stellar atmosphere (abundance stratification). 
These theoretical predictions can be directly tested through the 
comparison with empirical maps of chemical elements inferred 
from observations. For a small number of elements, includ- 
ing Ca, an effect of the vertically stratified element distribu- 
tion on the spectral line profiles was demonstrated in early stud- 
ies (Borsenberger at al. '1981). However, the absence of high- 
resolution, high signal-to-noise spectroscopic observations did 
not allow a robust comparison between observations and theo- 
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retical diffusion modelling. This step was carried out by Babel 
(119921 ). who calculated the Ca abundance distribution in the at- 
mosphere of magnetic Ap star 53 Cam and showed that the un- 
usual shape of Ca ii K line - a combination of the wide wings 
and extremely narrow core (Babel il994i Cowley et al. 12006!) - 
is a result of a step-like Ca distribution with abundance decrease 
at logT5ooo ~ -1- Following Babel, the step-function approxi- 
mation for the abundance distributions was employed in many 
stratification studies based on the observed profiles of spectral 
lines (Savanov et al.'200r, Bagnulo et al. lSOOTl Wade et al. l2003l 
Ryabchikova et al 2002 2005 2006; Glagolevskii et al. l2005l 
Cowley et al. 120071) . 

Ca was found to be stratified in nearly the same way as in 
53 Cam (enhanced concentration of Ca below logTsooo ~ -I 
and its depletion above this level) in all stars for which stratifi- 
cation analysis have been performed; /3 CrB (Wade et al. 20031, 
7 Equ (R yabchikova et al. l2002] l. HD 204411 (Ryabchikova et 
al.'2005'), HD 13379 2 (Koc hukhovet al. l2006al) and HD 144897 
(Ryabchikova et al. |2006l l. In addition to remarkable diffusion 
signatures, recently another Ca anomaly was detected, first in the 
spectra of the Hg-Mn stars by Castelh & Hubrig (12004 1) and then 
in Ap stars by Cowley & Hubrig ( 120051) . These authors found a 
displacement of the lines of Ca ii IR triplet due to significant con- 
tribution of the heavy Ca isotopes. Ryabchikova, Kochukhov & 
Bagnulo (see review paper by Ryabchikova (120051) ) were the first 
to apply spectrum synthesis calculations to investigate effects of 
Ca isotopes on the calcium line profile shape in magnetic CP 
stars. In particular, we have demonstrated the general anticor- 



relation between the presence of heavy Ca and magnetic field 
strength: in Ap stars the contribution of heavy Ca isotopes de- 
creases with the increase of the field modulus, and disappears 
when the field strength exceeds ~ 3 kG. Cowley et al. (2007) 
have come to a similar conclusion, but with reservations. 

In the present paper we summarise our detailed analysis of 
the vertical stratification of Ca abundance in the atmospheres of 
magnetic Ap stars of different temperatures and magnetic field 
strengths. We combine this vertical Ca mapping with the anal- 
ysis of the Ca isotopic anomaly and its dependence on height, 
based on the spectrum synthesis modelling of the IR triplet Ca ii 
line at A 8498 A. Our study is the first to present a homogeneous 
and systematic determination of the vertical stratification of a 
given element in the atmospheres of a large number of stars, en- 
abling us to study the signatures of atmospheric atomic diff'usion 
in the presence of strong magnetic field as a function of stellar 
parameters and magnetic field intensity. 

Our paper is organized as follows. In Sect. |2] we de- 
scribe spectroscopic observations and the data reduction. 
Determination of the stellar atmospheric parameters is detailed 
in Sect. [3] Sects. |4]-|6] present our methodological approach to 
magnetic spectrum synthesis, determination of Ca stratification 
and the study of Ca isotopic composition. Our findings are sum- 
marised in Sect. Hand are discussed in Sect.|8] 

2. Observations and data reduction 

Twenty-three slowly rotating Ap stars were chosen for the Ca 
stratification analysis. The list of program stars is given in 
Table [1] In addition, three stars, Procyon (=HD 61421, spec- 
ti-al type F5), HD 27411 (=HIP 20106, spectral type A3m), and 
HD 73666 (=40 Cnc, spectral type Al V) were used as standards 
for the Ca isotopic study. Of the 23 magnetic Ap stars included 
in our study, 12 stars show high-overtone acoustic p-mode pul- 
sations and thus belong to the group of rapidly oscillating Ap 
(roAp) stars. 

We note that metallic line stars show signatures of strong tur- 
bulence and convection in their atmospheres (Landstreet ll998l 
Kochukhov et al. 2006bl l. Theoretical diff'usion models for Am 
stars (e.g. Michaud et al. |2005 ) achieve reasonable success in 
reproducing the surface abundance patterns assuming full mix- 
ing of the outer envelope. Thus, despite non-solar abundances, 
strong mixing inhibits development of the vertical chemical gra- 
dients in the line-forming region and for this reason Am star 
HD 2741 1 can be considered as a comparison object in the con- 
text of the present investigation. 

For the majority of our targets high-resolution, high signal- 
to-noise-ratio spectra were obtained with the UVES instrument 
(Dekker et al. |2000] | at the ESO VLT in the context of pro- 
gram 68.D-0254. The observations were carried out using both 
available dichroic modes (standard settings 346+580 nm and 
437+860 nm). In both the blue arm and the red arm the slit width 
was set to 0.5", for a spectral resolution of about 80000. The 
slit was oriented along the parallactic angle, in order to mini- 
mize losses due to atmospheric dispersion. Almost the full wave- 
length interval from 3030 to 10400 A was observed except for a 
few gaps, the largest of which was at 5760-5835 A and 8550- 
8650 A. In addition, there are several small gaps, about 1 nm 
each, due to the lack of overlapping between the echelle orders 
in the 860U setting. 

Spectra of HD 24712 and HD 61421 (Procyon) were ac- 
quired with the UVES at VLT as part of the UVESpop projecfl 

' |http : //www ■ eso . org/uvespop7 1 



(Bagnulo et al. l2003al) . The UVES observations of the magnetic 
Ap star HD66318 employed in our paper were described by 
Bagnulo et al. (I2003bl l. All three stars were observed using the 
same instrumental settings as the Ap targets from our main data 
set. Spectrum of HD 73666 was obtained with the ESPaDOnS 
spectropolarimeter at the Canada-France-Hawaii Telescope and 
analysed for abundances by Fossati et al. (I2007I I. A reduced 
spectrum was kindly provided to us by the authors. 

The UVES spectra have been reduced with the automatic 
pipeline described in Ballester et al. (2000). For all settings, 
science frames are bias-subtracted and divided by the extracted 
flat-field, except for the 860 nm setting, where the 2-D (pixel- 
to-pixel) flat-fielding is used, in order to better correct for the 
fringing. Because of the high flux of the spectra, we used the 
UVES pipeline average extraction method. 

Due to the gaps in spectral coverage, only one line of the Ca ii 
IR triplet, A 8498 A, could be observed in all stars and is acces- 
sible for modelling. This line is the weakest among the three IR 
triplet transitions, which facilitates quantitative analysis, espe- 
cially in cooler stars where the triplet Ca ii lines become partic- 
ularly strong. 

All three lines of the Ca ii IR triplet overlap with the hydro- 
gen Paschen lines, which become very prominent in the spectra 
of hotter Ap stars (see Fig.[T]). The 8498 A line is located further 
away from the centre of the nearest Paschen line (PI 6) than the 
other two lines of the IR triplet. The hydrogen line blending for 
this Ca II transition is also weaker than for the other two Ca ii 
lines. Moreover, being the weakest among three IR triplet lines 
the 8498 A line is less dependent on the accuracy of Stark and 
Van der Waals broadening parameters. These circumstances sug- 
gest the Ca II 8498 A line is most suitable for detailed modelling. 

The studied Ca ii line lies rather close to the gap in the wave- 
length coverage of the UVES setting employed. An additional 
difficulty arises due to the broad overlapping absorption pro- 
duced by the Paschen lines. These properties of the data compli- 
cate continuum normalization in the region around Ca ii 8498 A. 
In fact, with the available data, we are generally unable to rectify 
the spectra in the usual manner Instead, for each star observa- 
tions were adjusted to match calculated Paschen line spectrum 
in a small region around the P16 line. This procedure was car- 
ried out in two steps. First, a high-order spline function was em- 
ployed to trace the wings of the hydrogen line, and observations 
were divided by this fit. Secondly, observed spectra were mul- 
tiplied by the theoretical hydrogen line spectrum calculated us- 
ing individual stellar atmospheric parameters and hydrogen line 
opacity described in Sect. |4] 

For most program stars spectral coverage of our UVES 
spectra precluded analysis of the two stronger lines of the IR 
triplet, Ca ii 8542 and 8662 A. However, for one of the targets, 
HD 24712, we have investigated these lines using complimen- 
tary observations obtained with the ESPaDOnS instrument at 
CFHT. The ESPaDOnS spectrum HD 247 1 2 is the average of 8 1 
time-resolved observations of this star analysed by Kochukhov 
& Wade ( I2007I I. We refer to this paper for the details of reduc- 
tion of these data. Unhke the UVES spectrum of HD 247 12, the 
ESPaDOnS data were collected close to the phase of magnetic 
maximum (JD 2453744.71-2453744.83). 

3. Fundamental stellar parameters 

Fundamental parameters of the program stars are given in 
Table [T] For most stars effective temperatures Teff and surface 
gravities log g were taken from the recent studies (references are 



Table 1. Fundamental parameters of the program stars. The columns give the target HD number, range of Julian dates of its obser- 
vation with UVES, effective temperature, surface gravity, projected rotational velocity, estimate of the mean magnetic field modulus 
and the ratio of radial to azimuthal field strength adopted in spectrum synthesis. Rapidly oscillating Ap (roAp) stars are identified 
by asterisk. 
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provided in the last column of Table [B- For HD 965, HD 47103, 
HD 118022 and HD 134214 atmospheric parameters were de- 
rived using the Stromgren photometric indices extracted from 
the catalogue of Hauck & Mermilliod ( 19981 ). We used cal- 
ibrations by Moon & Dworetsky (1985) and by Napiwotzki 
et al. (119931 ) as implemented in the TEMPLOGG code (Rogers 
119951 ). For HD 1 1 1 133 the effective temperature was taken from 
Kochukhov & Bagnulo (120061 ). while other parameters were 
derived in the present study. For HD 75445, HD 176232 and 
HD 203932 effective temperatures were further refined by fitting 
the Wa profile. This resulted in a small correction relative to the 
published Teff values. 

The mean surface magnetic fields (BJ were derived us- 
ing spectral lines which exhibit resolved and partially resolved 
Zeeman splitting patterns. In all stars rotational velocities were 
estimated by modelling profiles of the magnetically insensitive 
Fe I lines at AA 5434.5 and 5576. 1 A. Model atmospheres for all 
program stars were calculated with the ATLAS9 code (Kurucz 
119931 ). employing opacity distribution functions (ODFs) with an 
enhanced metallicity. For stars with especially strong magnetic 
fields we used ODFs with a non-zero pseudo-microturbulence 
velocity in order to simulate the modification of the line opac- 
ity due to magnetic intensification of spectral lines. All spec- 
trum synthesis for magnetic stars employed zero microturbu- 
lence. For comparison stars we have adopted the following val- 



ues: = 1-8 kms ' for Procyon (Allende Prieto et al. 120021 ). 
1.9 km s"' for HD 73666 (Fossati et al.|2007l) and 2.5 kms"' for 
HD 2471 1 (determined in this paper). 

4. Spectrum synthesis calculations 

Analysis of the Ca stratification and isotopic composition was 
based on detailed magnetic spectrum synthesis calculations with 
the SYNTHMAG code (Kochukhov I2007al ). This program allows 
to solve numerically polarized radiative transfer equation and 
calculate theoretical stellar spectra in four Stokes parameters 
for the prescribed model atmosphere and vertical distribution 
of any number of chemical elements. The SYNTHMAG calcula- 
tions performed in this paper are based on a simplified model 
of the stellar magnetic field topology, characterized by a sin- 
gle value of the field modulus and a homogeneous field distri- 
bution over the surface of the star. When rotation broadening 
is significant or the field is too weak, we assume purely radial 
field orientation. In other cases the effective field orientation is 
inferred from the profiles of magnetically sensitive lines and is 
parameterized with the ratio of the radial to azimuthal field com- 
ponents (see Table [T]). Although undoubtedly rather simplified, 
this magnetic model is quite successful in explaining the shapes 
of resolved Zeeman split lines in the Stokes / spectra of many 
Ap stars (Kochukhov et al. i2002l Nielsen & Wahlgren l2002l 




8500 8600 
Wavelength (A) 

Fig. 1. Synthetic spectra calculated with different treatment 
of the overlapping hydrogen line opacity. Usual calculation 
(HLINOP routine, dashed line) is compared with spectrum syn- 
thesis based on the occupation probability formalism (HBOP rou- 
tine, solid line). In both cases ATLAS9 model atmosphere with 
Teff - 8500 K, and log g - 4.0 was adopted. The vertical bars 
show position of the Can IR triplet lines. 



Leone et al. 120031 Ryabchikova et al. 120061 ) and, in some cases, 
performs even better than low-order multipolar field geometries 
(Kochukhov|2007b|l- 

In all spectrum synthesis calculations in this paper we use 
zero microturbulence because the SYNTHMAG code models mag- 
netic intensification by directly including magnetic field effects 
in the solution of the polarized radiative transfer equation. The 
true microturbulence is thought be absent in Ap stars because 
kG-strength magnetic field is quite efficient in suppressing con- 
vective motions. 

Here we analyse the target stars assuming a homogeneous 
distribution of Ca over their surfaces. This assumptions appears 
to be reasonable given a slow rotation of the majority of the tar- 
gets. Moreover, no evidence of substantial horizontal Ca abun- 
dance gradients exist for any of the studied stars, although many 
of them were observed at high resolution more than once (e.g., 
Mathys et al. 119971 Ryabchikova et al. I2004bl) . The only star 
from our sample for which the Ca surface inhomogeneity was 
investigated previously, is HD 24712 (Liiftingeret al. 2007 ). The 
full range of the horizontal Ca abundance variation in this star is 
only 0.2 dex. The resulting effect on Ca the line profiles is neg- 
ligible compared to the vertical abundance jump of more than 
3 dex (see below). We also note that, even if hypothetical high- 
contrast Ca spots exist on the surface of some of our stars, their 
effect will be to weight our modelling results to specific surface 
areas. But none of the potential horizontal inhomogeneity ef- 
fects can possibly mimic stratification or isotopic anomaly stud- 
ied here. 

An important modification to the SYNTHMAG treatment of the 
hydrogen line opacity was introduced to improve analysis of the 
IR Can triplet lines. These Can transitions blend with the high 
members of the Paschen series, which become prominent in Ap 
stars with T^ff > 8000 K. The usual calculation of the overlapping 
hydrogen line wings, which adds linearly opacities due to bound- 
free and bound-bound hydrogen opacity, becomes increasingly 
inaccurate as one approaches the series limit. As illustrated in 
Fig.[Tl this leads to a spurious reduction in calculated fluxes. The 
correct treatment of the overlapping hydrogen line and continu- 



ous opacity should be done according to the occupation proba- 
bility formalism (Daeppen et al. ll987i Hubeny et al. 1994). This 
technique was implemented in SYNTHMAG with the help of HBOP 
procedure, developed by R Barklem at Uppsala Observator}Q. 
This hydrogen opacity code, which also incorporates previous 
developments in the hydrogen line broadening theory (Stehle & 
Hutcheon ll999l Barklem et al. 120001 ). was applied for the spec- 
trum synthesis of all stars in our sample. Fig. [1] shows that cor- 
rect treatment of the hydrogen opacity significantly reduces pre- 
dicted Paschen line absorption at the position of the Ca ii 8498 A 
line analysed here. 

Magnetic spectrum synthesis with SYNTHMAG was cou- 
pled with the vertical abundance mapping procedure DDAFIT, 
described by Kochukhov (2007a) and used previously by 
Ryabchikova et al. (2005 2006 ). This routine provides a graphi- 
cal and optimization interface to SYNTHMAG, allowing the user to 
find parameters of a simple stratification model by fitting a large 
number of spectral lines of a given element. 

5. Ca stratification analysis 

Before performing detailed study of the IR Ca ii A 8498 A line 
profile, one has to derive vertical distribution of Ca abundance 
in the atmospheres of Ap stars. Without this step, a quantitative 
analysis of isotopic anomaly is practically impossible because 
profiles of all strong Ca lines, including the IR triplet, are dra- 
matically affected by the vertical stratification of this element. 
In all program stars Ca stratification was derived using a set of 
spectral lines in the optical region, for which no indication on the 
significant isotopic shifts exists. In this way we decouple strati- 
fication analysis from the study of isotopic shifts. At this stage 
we also avoid using the Ca ii 8498 A line because of the possible 
sensitivity to the errors in effective temperature due to the special 
normalization procedure adopted for this line (see Sect. 2). 

Atomic parameters of the Ca lines are given in Table|2] 
Stratification analysis requires high accuracy not only for the os- 
cillator strengths but also for the damping parameters, because 
Ca has a tendency to be concentrated close to the photospheric 
layers where the electron density is high. Accurate treatment of 
damping wings is particularly important for the lines of ionized 
Ca. For the Ca ii lines at AA 3158, 3933, 8248, 8254, 8498, 8912, 
8927 A the Stark damping constants were taken from the paper 
by Dimitrijevic & Sahal-Brechot ( 119931 ). where semi-classical 
calculations as well as a compilation of the experimental data 
was presented. For the rest of Ca lines the Stark damping con- 
stants calculated by Kurucz (119931 ) were adopted. The oscillator 
strengths were taken mostly from the laboratory experiments - 
these data were thoroughly verified by the recent NLTE analysis 
of calcium in late-type stars (Mashonkina et al. 2007 ). Because 
of the large range in effective temperatures and magnetic field 
strengths, we could not use the same set of lines for all stars. 

The Ca stratification analysis was performed using the step- 
function approximation of the abundance distribution (for details 
see, e.g., Ryabchikova et al. 2005). This parameterized model of 
the vertical stratification is characterized by the upper and lower 
abundance values, as well as by the position and width of abun- 
dance jump. In each star these parameters are constrained by si- 
multaneous fit to the profiles of many lines. These observational 
data are sufficient to constrain stratification parameters with a 
good precision. 

We start stratification analysis of Ca by computing the op- 
timal homogeneous Ca abundance for a chosen set of spectral 
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Fig. 2. Comparison between the observed line profiles (symbols) of HD 176232 (lOAql) and theoretical spectra calculated with 
the best-fitting stratified Ca abundance distribution (solid line) and with the homogeneous Ca abundance (dashed line). The plot- 
ted spectral windows show segments of the HD 176232 spectra used for determination of Ca stratification; symbols show actual 
observed points. 



Table 2. Atomic parameters of the spectral lines employed in the 
Ca stratification analysis. The columns give the ion designation, 
central wavelength of the transition, the excitation potential of 
the lower level, oscillator strength, the Stark damping constant, 
and the reference for the oscillator strength. 
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-0.079 


-3.70 


SMP 


Cai 


5857.451 


2.933 


0.240 


-5.42 


S 


Cai 


5867.562 


2.933 


-1.570 


-4.70 


S 


Cai 


6122.217 


1.896 


-0.316 


-5.32 


SO 


Cai 


6162.173 


1.899 


-0.090 


-5.32 


SO 


Cai 


6163.755 


2.521 


-1.286 


-5.00 


SR 


Cai 


6166.439 


2.521 


-1.142 


-5.00 


SR 


Cai 


6169.042 


2.253 


-0.797 


-5.00 


SR 


Cai 


6169.563 


2.256 


-0.478 


-4.99 


SR 


Cai 


6449.808 


2.521 


-0.502 


-6.07 


SR 


Cai 


6455.598 


2.523 


-1.340 


-6.07 


S 


Can 


6456.875 


8.438 


0.410 


-3.70 


SMP 


Cai 


6462.567 


2.523 


0.262 


-6.07 


SR 


Cai 


6471.662 


2.526 


-0.686 


-6.07 


SR 


Can 


8248.796 


7.515 


0.556 


-4.60 


SMP 


Can 


8254.721 


7.515 


-0.398 


-4.60 


SMP 


Can 


8498.023 


1.692 


-1.416 


-5.70 


T 


Can 


8912.068 


7.047 


0.637 


-5.10 


SMP 


Can 


8927.356 


7.050 


0.811 


-5.10 


SMP 



Rerefences. T - Teo dosiou ( fT989t ; SG - Smi th & Gallagher {T966> ; 
SMP - Seaton et al. (TW); S - Smith ( fT988] l: SO - Smith & O'Neil 
([1975); SR - Smith & Raggett U98ll . 



lines. Then we let the DDAFIT procedure to vary parameters of 
the step-function until the adequate fit to the observed line pro- 
files is achieved. The outcome of this procedure is illustrated 
in Fig. |2] where results of the stratification analysis for the set 
of optical lines in HD 176232 (10 Aql) are displayed. Figure [3] 



shows that stratified Ca distributions obtained for HD 176232 
and other cool Ap stars also dramatically improve the fit to 
the strongest Can 3933 A line compared to the spectrum syn- 
thesis with a homogeneous Ca distribution. In Figs. |2] and |3] 
synthetic profiles calculated with the uniform Ca distribution 
log(Cfl/Mot) - -5.14 are shown by dashed lines while those cal- 
culated with the best-fitting stratified Ca distribution are shown 
by the solid lines. Clearly, even a very schematic Ca strati- 
fication model used here provides a substantial improvement 
of the agreement between observations and theoretical calcula- 
tions. We use standard deviation to characterize the quality of 
the fit. This parameter is defined as a rms of the observed minus 
calculated spectrum, evaluated on the wavelength grid of obser- 
vations within the spectral windows covering lines of interest 
(see Fig. For HD 176232 and most other stars the stratified 
Ca abundance yields roughly two times smaller standard devi- 
ation compared to the homogeneous Ca distribution. The final 
Ca stratification model inferred for HD 176232 is illustrated in 
Fig.El 

To verify that the step-function approximation provides a 
realistic description of the Ca distribution we performed for a 
few stars with weak magnetic fields (HD 133792, HD 176232, 
HD 203932) stratification analysis using the Vertical Inverse 
Problem code, VIP. This program does not use any a priori 
assumptions about the shape of the vertical element abundance 
profile but currently cannot be applied to stars with a strong mag- 
netic field (Kochukhov et al. l2006al l. A comparison between the 
Ca distributions derived for HD 176232 with DDAFIT and with 
VIP is shown in Fig. |4] The same level of agreement was ob- 
tained for the two other stars. 

However, despite general success of the DDAFIT modelling, 
in a few stars the step-function approximation could not provide 
an adequate description of the full set of spectral lines. The ob- 
vious reason is the use of normal non-magnetic stellar model 
atmosphere with a homogeneous element distribution for a star 
with abundance stratification, and a deviation of the Ca abun- 
dance distribution from the simple step-function form. 
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Fig. 3. Comparison between the observed profiles (symbols) of 
the Can 3933 A line and calculations with the stratified (solid 
line) and homogeneous (dashed line) Ca distributions for 3 pro- 
gram stars. The spectra of HD 247 12 and HD 176232 are shifted 
downwards for display purpose. 



In cooler stars we do not systematically use the strongest 
Can 3933 A line in the stratification analysis due to difficulties 
of treating blends in the extended wings and establishing contin- 
uum. The range of the formation depths of the remaining weaker 
optical lines is different from the IR triplet lines of interest, there- 
fore Ca abundance in the upper atmospheric layers is poorly con- 
strained by the optical lines and may not be accurate enough for 
the description of the cores of the IR Can triplet lines. This is 
often reflected by a rather large formal error of the abundance in 
the upper layers compared with the corresponding error of the 
abundance in the lower layers. For several stars (see Table |4]i a 
change of the Ca abundance above in the upper layers within 3cr 
was introduced to fit overall intensity of the IR triplet line cores. 
In addition for a few other stars further improvement of the fit 
quality was possible by introducing a more complicated strati- 
fication profile starting from the results obtained with DDAFIT. 
This step was performed manually for 5 out of 23 Ap stars. In 
all the cases when we have modified Ca stratification determined 
with DDAFIT, the consistency of the altered vertical abundance 
profile with all Ca diagnostic lines was verified. 

The formal stratification analysis of the Am star HD 27411 
(whose Tgff is identical to that of HD 176232) and normal star 
HD 73666 (reff=9382 K) yields a 0.2 dex abundance jump. 
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Fig. 4. Vertical stratification and isotopic separation of Ca de- 
rived for HD 176232. The overall Ca abundance profile, inferred 
from the stratification analysis of optical lines, is shown with the 
solid line. The hatched areas demonstrate vertical separation of 
the light and heavy Ca isotopes required to fit the Ca n 8498 A 
line. The dashed line shows the Ca distribution derived the with 
VIP code using the same set of lines. 



which we consider to be insignificant. The average Ca abun- 
dance log(Ca/A^tot) = -5.63 in both stars represents very well 
most of the lines, including the resonance Can 3933 A line. 
In both stars 0.1 dex lower abundance is required to fit the 
Can 3933 A profile compared to Can 8498 A line. This inter- 
nal agreement is excellent considering difficulties in continuum 
normalization. It confirms that the Ca abundance variation by 
several dex over the vertical span of the Ap-star atmospheres is 
real. 



6. The Ca isotopic anomaly 

Ca has six stable isotopes with atomic numbers 40, 42, 43, 44, 
46, 48. In the terrestrial matter Ca mixture consists mainly of 
"^''Ca (96.9 % according Anders & Grevesse'T989). Table[3]gives 
wavelengths of all Ca isotopes in the Can 8498 A transition 
based on the isotopic shifts measured by Nortershauser et al. 
(il998i ). We also list fractional oscillator strengths corresponding 
to the terrestrial isotopic mixture. 

Isotopic shifts were also measured for the Can resonance 
fines at AA 3933 and 3968 A (Martensson-Pendrifl et al. [T992] l. 
The wavelength shift between the ^^Ca and ^*'Ca isotopes is 
-0.009 A, and can be observed only for the Can 3933 A line 
core in hot stars with an especially large contribution of the "*^Ca 
isotope. One star from our sample, HD 133792, shows an iso- 
topic shift in the Can 3933 A line. 

Using the terrestrial isotopic mixture, we calculated the pro- 
file of the Ca n 8498 line in the spectra of our reference stars, 
Procyon, HD 2741 1 and HD 73666. Comparison of the theoret- 
ical spectra and observations is illustrated in Fig. |5] Although 
in the Procyon spectrum our LTE calculations cannot provide a 
very good fit to the deepest part of the core, no wavelength shifts 
are detected in either star. At the same time, the observed pro- 
file of the Ca n line in the spectrum of one of our program stars, 
HD 217522, presented in Fig.|5]exhibits a complex structure and 
is clearly redshifted, with the strongest component coinciding in 



Table 3. Atomic data for the isotopic components of Can 
A 8498 A. The fractional abundance of Ca isotopes, s, corre- 
sponds to the composition of the terrestrial matter 



A (A) 


isotope 


log gfs 


8498.023 


40 


-1.43 


8498.079 


42 


-3.60 


8498.106 


43 


-4.29 


8498.131 


44 


-3.10 


8498.179 


46 


-5.81 


8498.223 


48 


-4.14 
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Fig. 5. Comparison between the observed line profiles (symbols) 
of the Ca n 8498 A line and calculations for the normal Ca iso- 
topic mixture (solid line) in the spectra of Procyon, Am star 
HD 27411, and Al V star HD 73666. The observed spectrum 
of the cool (HD 217522) and tepid (HD 133792) Ap stars, pre- 
sented in the middle and bottom, clearly shows a displaced Ca ii 
line core. The spectra are shifted downwards for display purpose. 



wavelength with the expected position of the heaviest Ca iso- 
tope. This is a signature of the Ca isotopic anomaly and, as we 
will show below, of the vertical separation of Ca isotopes. 

The core of the IR Ca ii 8498 A line is formed higher than 
any of the optical lines, except Can 3933 A. For most stars the 
Can 3933 A line was not modelled in our stratification calcula- 
tions. Therefore, Ca abundance in the upper atmosphere of some 
stars may be somewhat uncertain, as all other optical lines are 
not sensitive to the abundance variations above log rsooo - -2.0 
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Fig. 6. Derivation of the isotopic composition and its height de- 
pendence for the Can 8498 A line in HD 176232 (lOAql). The 
three panels show the UVES observations (full circles) com- 
pared to various modelling attempts, as follows, a) Dashed line: 
theoretical spectrum obtained with a homogeneous Ca distribu- 
tion. Solid line: theoretical spectrum calculated for a stratified 
Ca abundance. A terrestrial mixture of Ca isotopes is adopted 
for these computations, b) Solid line: theoretical spectrum ob- 
tained for a terrestrial Ca isotopic composition (same as panel 
a); Dashed line: theoretical spectrum obtained for a 50:50 ratio 
of "^''Ca to '*^Ca-i-'***Ca. c) Solid Une: the Ca ii line profile com- 
puted for the best-fitting model of vertical isotopic separation 
(see Fig.HJl. Dashed and dash-dotted lines show contributions of 
''"Ca and "^*Ca-i-^^Ca, respectively. The wavelength position of 
^°Ca and "^^Ca is indicated in each panel. These plots show that 
an acceptable fit to the observations can only be achieved with a 
depth-dependent Ca isotopic mixture. 



to -2.5. However, for the majority of our targets the Ca abun- 
dance in the upper atmosphere, responsible the strength of the 
Ca II 8498 A Hne core, is defined by the slope of the abundance 
gradient in the jump region. If the Ap atmospheric structure is 
close to the normal ATLAS9 atmosphere adopted in our analy- 
sis, then the Ca ii 8498 A line should be well-described by the 
Ca abundance distribution derived from optical lines. Our cal- 
culations show that while it is correct for the observed total in- 
tensity, in many of the program stars we cannot fit the line cores, 
which are often redshifted. This situation is illustrated in Fig.|6^, 
where we compare the observed spectrum of HD 176232 with 
the synthetic spectrum calculated for the terrestrial Ca isotopic 
mixture and Ca abundance distribution (Fig.|4]i derived from op- 
tical lines. One immediately notices that while the line wings are 
explained by our calculations, the line core cannot be fitted with 
the terrestrial Ca isotopic mixture. 

Cowley & Hubrig ( 20051 1 showed that the red-shift of the Ca 
IR line core arises due to the heavy Ca isotopes and claimed that 
"in extreme cases the dominant isotope is the exotic '**'Ca". A 
simple interpretation of the anomaly observed in the Ca ii 8498 A 
line core is to suggest that heavy Ca isotopes are strongly en- 
hanced and even dominant throughout the atmospheres of some 
magnetic Ap stars. However, our magnetic spectrum synthesis 
calculations demonstrate that this hypothesis is incorrect. On the 
example of HD 176232 we can show that using a model with 
stratified Ca distribution found previously and assuming a depth- 
independent enhancement of "^^Ca and "^^^Ca relative to '^'^Ca, the 
center-of-gravity of the Ca ii 8498 A line core is shifted to the 
red, but the overall fit to the line profile does not improve con- 
siderably (Fig. I6J5). In disagreement with observations, a con- 
stant overabundance of the heavy Ca isotopes gives rise to the 
red-shifted long-wavelength wing of the Ca 11 line. 

Here we propose a different explanation of the Ca 11 8498 A 
line shape. Observations of HD 176232 and other stars with dis- 
placed Ca line core show a shallow line with wide wings at the 
position of '^^Ca isotope, which means that this line component 
is formed in the lower atmospheric layers. On the other hand, at 
the position of '*^Ca we see a sharp deep line, which is responsi- 
ble for a characteristic steep intensity gradient in the red wing of 
the Ca II 8498 A line. We observe the same picture in all cool Ap 
stars with significant isotopic shifts of the Ca 11 line core. Given 
the Ca stratification inferred for all Ap stars, this deep compo- 
nent lacking developed wings can be explained by the absorption 
in the upper atmosphere, above the Ca abundance jump or in the 
outermost part of the transition region. Thus, observations of the 
Can 8498 A line in Ap stars should be interpreted in terms of 
vertical separation of Ca isotopes in the stellar atmosphere. 

If we separate ""^Ca and the sum of '*''Ca and ^^Ca isotopes in 
the atmosphere of HD 176232 as indicated in Fig.|4](the bound- 
ary between heavy and light Ca is at log 75000 = - 1 -2), a satis- 
factory agreement between the observed and calculated spectra 
is obtained (see Fig.|6]:). The dominant atmospheric constituent 
is still the normal isotope '^''Ca, but the shape of the Ca 11 8498 A 
line core is influenced by the high-lying cloud of heavy Ca, in 
which "^^Ca and '^'^Ca contribute 35% and 65%, respectively. 

Thus, in our modelling of the Ca isotopic anomaly in Ap 
stars we adopt an initial model in which the total element abun- 
dance is given by the stratification profile obtained previously, 
while the vertical separation of the light and heavy Ca isotopes 
occurs nearly instantaneously and is characterized by a single 
transition depth, deduced from the Can 8498 A line. Relative 
contributions of '^''Ca and '^^Ca to the heavy Ca layer in the up- 



per atmosphere are also adjusted based on the Can IR triplet 
line. 
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Fig. 7. Comparison of the observed IR triplet lines (symbols) in 
the ESPaDOnS spectrum of HD 24712 and theoretical spectrum 
synthesis for the best-fitting stratified Ca abundance distribution 
(solid line) with Ca isotopic separation and with the stratified 
distribution of "'^^Ca only (dashed line). 

Parameterization of the vertical separation of Ca isotopes 
that we adopt here is just one of several models possible. For 
instance, we find that an acceptable fit to the Can 8498 A line in 
some of the target stars could be obtained with a model where 
"'^^Ca is homogeneously distributed while ''"Ca is stratified. (Such 
a model is achieved by placing the border between Ca isotopes 
horizontally rather than vertically in Fig. HI Note, that formally 
'^'^Ca is not anomalous: its relative-to-hydrogen abundance in 
all stars does not exceed an abundance in the normal solar Ca 
mixture.) This model turns out to be nearly equivalent to the 
one we use because the boundary between the heavy and light 
Ca layers is usually located very close to the position of the 
Ca abundance jump. The "^''Ca contribution still dominates in 
the lower atmosphere, whereas ^'^Ca contributes to the light ab- 
sorption in the upper layers. However, for several stars where 
we need Ca enhancement in the uppermost layers, the model 
with homogeneously distributed heavy Ca isotopes clearly yields 
worse agreement with observations in comparison with the ver- 
tical separation model adopted in here. 

Our isotopic stratification model is based upon only one of 
the IR triplet lines. High-quality observations of the two other 
IR triplet lines are unavailable for the majority of our targets. 
For one of the cool Ap stars, HD 24712, we assessed the qual- 
ity of the fit to all three Ca ii lines using the spectrum obtained 
with ESPaDOnS at CFHT These data were collected at the ro- 
tational phase diff'erent from that of UVES observations, there- 
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Fig. 8. Ca stratification in the atmospheres of selected program 
stars. Upper panel: distribution of Ca in HD 66318, HD 75445, 
HD 134214, HD 166473 and HD 111133. For stars enriched in 
heavy Ca the filled circle indicates the bottom of the heavy Ca 
isotope layer. Lower panel: same for stars with complex dis- 
tribution of Ca abundance with height (HD965, HD 116114, 
HD 137909, HD 137949, HD 217522). 



fore we have performed an independent Ca stratification and iso- 
topic analysis for HD 24712 based on the ESPaDOnS spectrum 
alone. Fig.|2]compares calculated and observed Can 8498 A and 
8662 A line profiles. Due to specific position of the Ca ii 8542 A 
line near the beginning of the spectral order, which makes con- 
tinuum normalization uncertain, we did not include this line to 
the plot. 

Of course, the proposed model of Ca isotopic separation is 
at best a crude approximation of potentially much more complex 
isotopic stratifications. However, we find that, for many stars, 
introducing a combination of stratification and isotopic segrega- 
tion noticeably improves the fit to the IR triplet lines and thus 
provides a direct evidence for the Ca isotopic separation in the 
atmospheres of Ap stars. 

Finally, we note that the continuum normahzation in the re- 
gion around IRT lines is indirectly based on the adopted effective 
temperatures, which may introduce an additional uncertainty and 
sometimes lead to a poor fit in the line wings. However, we have 
verified that even considerable errors in Teff lead only to a small 
change of the overall depth of the Ca ii 8498 A line computed 
for a given Ca stratification. Since this line is located far enough 
from the H line core, its profile, especially the shape of the inner 
core region used in the isotopic analysis, is generally not aff'ected 
by the T^ff errors. 



7. Results 

7.1. Ca stratification 

The Ca stratification analysis described in Sect. |5] was applied 
to all Ap stars included in our sample. The presence of a strat- 
ified Ca distribution is inferred for all stars except the hot Ap 
star HD 170973, where stratification appears to be negligible. 
Table |4] gives parameters of the step-function model distribu- 
tions of Ca abundance. These vertical abundance profiles are 
characterized by an abundance jump in the atmospheric region 
-1.3 < logTsooo ^ -0.5, a 1-1.5 dex overabundance deep in the 
atmosphere and a strong Ca depletion above logrsooo ~ -1.5. 
The formal uncertainty of the stratification parameters is typi- 
cally 0.3-0.7 dex, 0. 1-0.3 dex and < 0. 1 dex for the upper abun- 
dance, lower abundance and the position of the Ca abundance 
jump, respectively. The Ca stratification in the atmospheres of 
selected program stars is shown in Fig. [8] 

For several stars a more complex vertical distributions of 
Ca abundance had to be introduced. Observations of HD 965, 
HD 137909, HD 137949 are better reproduced with an increase 
of Ca concentration in the upper atmospheric layers above 
log 

''"5000 — ~5. Although seemingly complicated, this picture 
does not contradict theoretical Ca diffusion calculations. Both 
Borsenberger et al. (1981. Fig. 6) and Babel (1992) obtained 
Ca abundance increase in the uppermost layers, above the main 
abundance jump. However, NLTE treatment of the Ca lines 
formation is needed to investigate stratification of Ca in the 
upper atmospheric layers in more detail. For two other stars, 
HD 116114 and HD 217522, a structure of the transition zone 
somewhat more complex than a linear transition between two 
constant values is required to fit observations. Figure [8] (lower 
panel) illustrates the final chemical stratification profiles adopted 
for the five stars where Ca stratification deviates from a step- 
function profile. 

We note that for all stars the DDAFIT inversion was per- 
formed on the original column mass (logpx) depth grids of the 
ATLAS9 models. Table H gives vertical stratification parameters 
converted to the log tsooo scale. In the cases when transition zone 
is wide, extending to the lower atmospheric layers, transforma- 
tion between the column mass and the standard optical depth 
scale becomes non-linear, which leads to a distortion of the step- 
function shape (e.g., Ca distribution for HD 137949 in Fig.[8]l. 
For the same reason in some of the stars where this problem 
was encountered (e.g., HD 137949, HD 188041) the formal so- 
lution for the Ca distribution yields an unrealistically high Ca 
abundance in the lower atmosphere. However, we find that in 
these objects the actual inferred stratification profile often has a 

plateau with log(Cfl/Mot) ~ -3.5 in the layers around log T5000 ~ 
0. This value of the abundance gives the actual element concen- 
tration that influences the Ca line formation. The formal high Ca 
abundance occurs far below the photosphere (e.g., the bottom 
panel of Fig.[8]shows this situation for HD 137949). It does not 
influence the observed Ca line profiles. Therefore, for all stars 
marked by a letter 'c' in Table |4] we took the Ca abundance at 
log T5000 ~ as log(Cfl/Mot)io, modifying the position and the 
width of the abundance jump correspondingly. 

Homogeneous determination of the calcium vertical distribu- 
tion in a large number of magnetic Ap stars allows us to search 
for possible dependence of diffusion signatures on the stellar at- 
mospheric parameters and magnetic field strength. However, ex- 
cept of the marginal tendency of the Ca abundance in the up- 
per layers to be higher in hotter stars, we found no clear cor- 
relations between each of the four stratification parameters on 
the one hand and Teff and (B^) on the other hand. Although the 



Table 4. Parameters of the Ca stratification and isotopic anomaly determined for magnetic CP stars. The columns give HD number 
of stars, step-function parameters of the vertical Ca distributions (abundance in upper atmosphere, abundance in the lower layers, 
position of the abundance jump and the width of the transition zone), parameters of the Ca isotopic separation if applicable (the 
optical depth of the boundary between the heavy and light Ca layers, the dominant isotopic component of the heavy Ca cloud). The 
last column indicates deviation of the final Ca stratification profile from a simple step-function. 



HD 




Stratification parameters 




Ca isotopic anomaly 


Comment 


number 


log(Ca/A'tot)up 






A log Tsoooistep) 




isotope 




965 


-8.4 


-4.9 


-0.7 


0.1 


-4.9 


■*'^Ca 


a, d 


24712 


-8.7 


-5.1 


-1.2 


2.5 


-1.4 


46ca 


c,d 


29578 


-8.9 


-4.3 


-1.2 


0.7 








47103 


-7.7 


-4.6 


-0.8 


0.3 








66318 


-8.0 


-5.0 


-0.6 


1.1 








75445 


-8.9 


-4.9 


-0.9 


0.2 


-2.0 


"^'Ca 




101065 


-8.9 


-5.3 


-0.4 


0.4 


-0.2 


"^Ca 




111133 


-10.0 


-6.3 


-2.0 


0.4 


-1.5 


"^Ca 




116114 


-7.4 


-4.7 


-1.3 


0.2 


-3.8 


"^^Ca 


b 


118022 


-7.0 


-2.0 


-0.1 


0.3 


-1.5 


"'^Ca 




122970 


-7.7 


-5.1 


-1.5 


2.8 


-1.6 


"^Ca 


c,d 


128898 


-8.5 


-4.0 


-1.0 


2.0 


-1.9 


48 Ca 


c,d 


133792 


-8.1 


-5.6 


-0.6 


0.1 


-1.0 


48 Ca 




134214 


-8.4 


-4.8 


-0.9 


0.9 


-1.5 


48 Ca 




137909 


-8.2 


-4.4 


-1.0 


0.4 


-4.4 


«Ca 


a 


137949 


-9.6 


-3.5 


-1.0 


1.8 


-4.9 


«Ca 


a, c 


144897 


-8.5 


-5.2 


-1.9 


0.5 








166473 


-8.5 


-3.8 


-1.0 


1.1 






d 


170973 


-5.3 


-5.0 


-1.1 


0.1 








176232 


-9.0 


-4.1 


-1.2 


1.8 


-1.2 


''^Ca 


d 


188041 


-7.4 


-3.1 


-1.2 


2.5 


-3.7 


**Ca 


c\d 


203932 


-8.7 


-4.5 


-1.2 


2.6 


-2.1 


«Ca 


c,d 


217522 


-8.8 


-4.8 


-0.3 


0.7 


-0.6 


''^Ca 


b 



Comments, a - abundance increase in the upper layers; b ~ complex transition zone; c - distortion of the step-function after logpx to logTsooo 
transformation; d -a. change in the upper abundance value. 



qualitative form of the Ca stratification is the same in all stars, 
possibly indicating the action of a universal radiative diff'usion 
process, the observed star-to-star variation in stratification pro- 
files appears to be irregular and unrelated to any other stellar 
characteristic. 

7.2. The Ca isotopic anomaly 

The Ca isotopic analysis procedure outlined in Sect. |6] was ap- 
plied to all stars included in our paper. Table |4] lists the lower 
boundary (uncertain to within ^0.1 dex) of the heavy Ca layer in 
the atmospheres of magnetic Ap stars for which the Ca ii 8498 A 
line could not be reproduced with the terrestrial isotopic mixture. 
Our estimate of the dominant contributor among different heavy 
Ca isotopes is also given. Out of 23 program stars, the presence 
of heavy isotopes is established in 17 objects. 

For three stars, HD965, HD 137909, HD 137949, the mod- 
elling of the IR triplet line requkes an increase of the Ca abun- 
dance in the upper atmosphere, above the abundance jump given 
by the simple stratification analysis. In these cases we attribute 
the high-lying cloud to heavy Ca. 

Figure |9] compares theoretical spectrum synthesis calcula- 
tions with observations for a subset of 12 cooler (T^ff < 9000 K) 
stars with different effective temperatures and different magnetic 
field intensity. The stars are arranged in order of increasing mag- 
netic field strength. It is evident that the presence of heavy Ca 
isotopes inversely coiTelates with magnetic field strength. This 
effect is further illustrated in Fig.fTTh. where we quantify the de- 
tection of heavy Ca by the ratio of calculated equivalent widths 
of the '*'*CaH-"^^CaH-"^**Ca to ""'Ca line components. In the stars 



with small to moderate magnetic fields we clearly see a signif- 
icant contribution of the heavy isotopes ^''Ca and "^^^Ca (large 
equivalent width ratio). This contribution decreases with the in- 
crease of magnetic field strength. Even in HD 137909 (fi CrB) 
with the mean magnetic field modulus {Bs)-5.4- kG one still 
needs a small contribution of "^^^Ca, but under the assumption of 
very specific Ca distribution shown in Fig. [8] No evidence for the 
presence of heavy Ca is found in cool Ap stars with (BJ ^6 kG. 

Due to the temperature dependence of the Ca line intensity, 
we used a smaller number of lines for stratification analysis of 
hot Ap stars (Tgff > 9000 K). In these stars stratification is de- 
fined mainly by the Ca ii lines. Among the neutral Ca lines only 
the resonance Ca 1 4227 A line could be utilized. The agreement 
of calculated and observed line profiles is poorer for hot stars. 
However, the general anticorrelation between the presence of the 
heavy Ca isotopes and magnetic field strength is still apparent 
(see Fig. [TOl i. As in cool stars, no evidence for the presence of 
heavy Ca is found in hot Ap stars with (BJ ~6 kG. The only ex- 
ception is HD 170973, which has no observable magnetic field 
and shows no contribution of the heavy Ca isotopes. We note 
that HD 170973 is also the only star for which we find no evi- 
dence of Ca stratification. This star also possesses the highest Ca 
abundance among the hot Ap stars in our sample. 

Thus, it appears that the mechanism responsible for the accu- 
mulation of heavy Ca isotopes in the line-forming atmospheric 
region works more efficiently in stars with a small magnetic field 
strength. Figure [TT] illustrates several other interesting relations 
between the Ca isotopic anomaly and stellar parameters that we 
have revealed in this study. The optical depth of the layer sep- 
arating the heavy and light Ca clouds depends on the magnetic 
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Fig. 9. Comparison of the observed (full circles) and calculated (lines) profiles of the Ca ii 8498 A line in a subset of program stars. 
Solid line shows theoretical calculations for the terrestrial Ca isotopic mixture and the Ca stratification derived from optical lines. 
The dashed line presents spectrum synthesis for the models with vertical separation of the Ca isotopes. The stars are arranged in 
order of increasing magnetic field strength. 



field strength. In stars with more intense fields heavy Ca iso- 
topes tend to accumulate much higher in the stellar atmospheres 
(Fig.fTTb). There is also an anticorrelation of the heavy Ca con- 
tribution and the total equivalent width of the Ca ii 8498 A tran- 
sition (Fig.fTTl:). which implies that the efficiency of the isotope 
separation process depends on the total line intensity in cases 
where such separation exists. 

We find no difference between roAp and non-pulsating Ap 
stars (shown with different symbols in Fig. [TTl ) regarding the 
presence of heavy Ca excess in their atmospheres and its cor- 
relation with the magnetic field strength and line intensity. 

8. Discussion and conclusions 

In this paper we have investigated the vertical stratification and 
isotopic anomaly of Ca in a large sample of cool magnetic Ap 
stars. Our study is the first to address the interesting problem of 
the presence of heavy Ca isotopes in chemically peculiar stars 
with detailed polarized radiative transfer calculations, which 
take into account the effects of the magnetic field and chemi- 
cal separation in stellar atmospheres. We derive stratification of 
Ca for 23 Ap stars using a sample of Ca i and Ca ii lines dis- 
tributed over a broad spectral range. All but one program stars 
clearly show signatures of the Ca stratification, whereas compar- 
ison stars reveal no inhomogeneities in the vertical Ca distribu- 
tion when analysed with the same techniques and atomic data. 
Although our stratification modelling was based on a simplified 
step-function approximation and adopted a simplified homoge- 



neous model for the magnetic field geometry, the inferred param- 
eters of the vertical Ca distributions impose important observa- 
tional constraints for theoretical models of the radiative diffusion 
processes in stellar atmospheres. 

Analysis of the IR Ca triplet line Ca ii 8498 A provided in- 
formation on the relative contribution of different Ca isotopes. 
Spectrum synthesis calculations show that significant isotopic 
shifts observed in the core of Ca ii 8498 A cannot be attributed 
to the overabundance of heavy Ca isotopes throughout the whole 
stellar atmosphere. Instead, we show that the Ca line profile 
shape is consistent with the vertical separation of different Ca 
isotopes, with heavy Ca located in a cloud above the most abun- 
dant isotope ^°Ca. Even though the presence of heavy Ca is 
prominent in the line core, the normal Ca isotope dominates the 
line wings and is more abundant in the lower atmosphere where 
the total Ca abundance is also much larger Thus, our tentative 
model calls for a less extreme heavy Ca enrichment than was 
suspected in previous investigations limited to the centroid mea- 
surements of the Ca ii 8498 A line core. 

Recently different stratification of He isotopes was found by 
Bohlender (2005) in the analysis of ^He and related stars. That 
paper used a similar approach to modelling chemical stratifica- 
tion, but has only addressed He vertical distribution in hot non- 
magnetic chemically peculiar stars. 

We have successfully used the model with Ca stratification 
and vertical isotopic separation to explain the appearance of the 
Can 8498 A in all stars showing excess of heavy Ca isotopes. 
The prominent anticorrelation between the presence of heavy 
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Fig. 11. Results of the analysis of Ca isotopic composition in Ap stars, a) ratio of calculated equivalent width of the heavy and light 
Ca components of the 8498 A line vs. magnetic field strength, b) Position of the boundary between the light and heavy Ca layers 
vs. magnetic field strength, c) The equivalent width ratio vs. the total equivalent width of the Can 8498 A line, d) The equivalent 
width ratio vs. the optical depth of the boundary between the light and heavy Ca layers. In all plots different symbols correspond to 
Ap stars with detected presence of heavy Ca (filled symbols) and those showing no heavy Ca isotopes (open symbols), as well as to 
roAp (circles) and non-pulsating Ap (triangles) stars. 



Ca and magnetic field strength, first reported by Ryabchikova 
(1200 5 ). is confirmed and strengthened in the present paper We 
find that only stars with sufficiently weak field show traces of 
heavy Ca. According to our knowledge, this interesting relation 
is the only case when definite dependence of the chemical abun- 
dance characteristic on the magnetic field strength is found for 
Ap stars. Furthermore, we find that in stars with stronger fields 
the heavy Ca isotopes tend to accumulate higher in the atmo- 
sphere. 

According to the results of our study, pulsating (roAp) and 
non-pulsating Ap stars are not distinguished by the character- 
istics of their Ca isotopic anomaly. Both groups of stars are 
equally likely to show an excess of heavy Ca and follow the same 
trend of the heavy isotope contribution versus the magnetic field 
strength and the Can 8498 A line intensity. 

If the overall distribution of Ca abundance in the atmo- 
spheres of Ap stars follows the predictions of the radiatively 
driven diffusion, our results on the isotopic separation favour the 
light-induced drift (LID) as the main process responsible for this 
separation. According to Atutov & Shalagin ( 1988), LID arises 
when the radiation field is anisotropic inside the line profile. 
Such an anisotropy takes place for a line of the trace isotope, for 
instance '^^Ca and "^^Ca in the terrestrial calcium mixture, which 
is located in the wing of a strong line of the main isotope '^''Ca. 
The main isotope induces the drift velocity for other isotopes. If 
the trace isotope's line is located in the red wing of the line due 



to main isotope, the drift velocity is directed towards the upper 
atmosphere and the trace isotopes are pushed upwards. This is 
in agreement with observed vertical distribution of Ca isotopes. 
The Zeeman splitting changes the line shape and decreases the 
flux anisotropy for the line of trace isotope. When magnetic field 
becomes strong enough, ~5-6 kG, the flux anisotropy disappears 
and the isotopic separation is ceasing. Therefore, the observed 
Ca isotopic anomaly in magnetic stars may be qualitatively ex- 
plained by the combined action of the radiatively-driven diffu- 
sion and the light-induced drift. Theoretic study of a combina- 
tion of the radiatively-driven diffusion and LID was presented by 
LeBlanc & Michaud (1993 ) for He. These authors showed that 
LID accelerates separation of ^He from "^He in hotter CP stars. 
Detailed theoretical chemical diffusion calculations (LeBlanc & 
Moninj2004) should incorporate LID in order to test our hypoth- 
esis that this effect may be important for the chemical transport 
processes in cool Ap-star atmospheres. 
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